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We study two dimensional stripe forming systems with competing repulsive interactions decaying as r~ a . 

We derive an effective Hamiltonian with a short range part and a generalized dipolar interaction which depends 
on the exponent a. An approximate map of this model to a known XY model with dipolar interactions allows us 
to conclude that, for a < 2 long range orientational order of stripes can exist in two dimensions, and establish 
the universality class of the models. When a > 2 no long-range order is possible, but a phase transition in the 
KT universality class is still present. These two different critical scenarios should be observed in experimentally 
relevant two dimensional systems like electronic liquids (a = 1) and dipolar magnetic films (a = 3). Results 
from Langevin simulations of Coulomb and dipolar systems give support to the theoretical results. 
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Two dimensional isotropic systems in which a short range 
attractive interaction competes with a repulsive interaction de¬ 
caying as a power law of the form r~ a have been widely stud¬ 
ied [1-8]. These include, as physically relevant examples, 
the dipolar (r~ 3 ) and the Coulomb (r _1 ) interaction as the 
repulsive part of the total energy of the system. Dipolar in¬ 
teractions competing with exchange and uniaxial anisotropy 
arise, e.g. in ultra-thin ferromagnetic films with perpendicular 
anisotropy [9-11] while long range Coulomb interactions ap¬ 
pear in low dimensional electron systems and may be relevant 
to understand the low temperature phase behavior of doped 
Mott insulators, two dimensional quantum Hall systems and 
high T c superconductors [12-15]. It is well known that un¬ 
der certain conditions of relative strength of interactions and 
external parameters these systems develop modulated stripe¬ 
like structures in two dimensions which break space rota¬ 
tional symmetry, similar to classical liquid-crystal systems, 
giving rise to smectic, nematic and hexatic phases [16-18]. 
This analogy, based on the 180° rotational symmetry of stripe 
structures and elongated liquid-crystal molecules, allowed to 
apply well known results for liquid-crystal systems [19, 20] 
to predict the qualitative, and to some extent also quantitative 
phase behavior of many systems with modulated order param¬ 
eters. Nevertheless, when it is important to understand the 
true nature of the thermodynamic phases, the analogy between 
stripe forming systems and classical liquid-crystals should not 
be taken at face value. The basic units in liquid-crystals are 
elongated molecules. A given molecule typically interacts 
with its near neighbors and due to its elongated form a ro¬ 
tation in 180° of a single molecule does not alter the energy 
of the system. On the other hand, the smallest relevant scale 
of a stripe system is the modulation length. At this scale, a ba¬ 
sic cell can be considered as containing a single interface and 
then it is a dipole of opposite densities with an average linear 
size equal to the modulation length. It is important to note 
that such dipoles will not be, in general, elementary electric 


or magnetic dipoles, their character will depend on the nature 
of the density order parameter under consideration. Having 
clarified this point, in principle all realistic low energy con¬ 
figurations of the system can be built from these dipole cells. 
Clearly, a 180° rotation of a dipole does change the energy 
of the system and then cannot be considered a local symme¬ 
try. The system is only symmetric under global rotations of 
180°. Furthermore, when long range interactions are present, 
it is well known that the behavior of the systems may be very 
different from those with only short range interactions, which 
represent the vast majority of classical liquid-crystal systems. 
A study of the nature of low temperature phases of stripe 
forming systems should take these elements into account. 

Consider a coarse grain Hamiltonian in two dimensions of 
the form 

'H.[<!>(x)\ = ^ j d 2 x (V0(£)) 

+ d 2 x j d 2 x' <f){x) J(x — x')<j>{x') 

+ jjj f d 2 x V(<j>(x)), (1) 

where /3 = 1/ksT and V(tj>) = —1</> 2 + is a local poten¬ 
tial that could be seen as an entropic contribution and which 
exact form is not important to our work. The long range re¬ 
pulsive interaction has the form J(x) = J/\x\ a which allows 
to analyze in a unified way short range (large a ) and long 
range (small a) interactions. Physically relevant examples 
are the Coulomb interaction (a = 1) and the dipolar inter¬ 
action between out-of-plane magnetic moments (a = 3). It 
is well known that at low temperatures this kind of systems 
display stripe-like patterns in the form of spatial modulations 
of the density (f>{x) [1, 21-23] in a direction represented by 
a wave vector /c 0 - Low energy excitations of the stripes can 
be described in terms of a displacement field u(x) in the form 
<f>{x) = J2 n Qn cos(nkox + nkou(x)), where x is the average 
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direction of the modulation and kg stands for the modulus of 
kg. If u(x) varies smoothly in space it is possible to define a 
local wave vector kg + k 0 Wu(x). 

The effective Hamiltonian (1) when expressed in terms of 
u(x) has local and non-local parts T-L = 'Hi -KH„z (see Supple¬ 
mental Material). Expanding the local component to quadratic 
order in the fluctuation field u, it can be written in Fourier 
space as [19, 22, 24]: 

Ui = Hoi + ^ J -^^{ixkl + ^ v ky)u(k)u{-k) ( 2 ) 

where ~/ x and ~/ y are elastic coefficients which are simply re¬ 
lated to the parameters of the original Hamiltonian and 'Hoi 
represents the local contribution to the energy for an unper¬ 
turbed stripe. It is well known that this form for the local 
fluctuations of the stripe pattern leads to a divergence of the 
mean square of the displacement field, implying the abscence 
of long range positional order in the system. This is the stan¬ 
dard situation in liquid-crystalline systems. We go on to con¬ 
sider the effect of the tail of the long range interaction in the 
fluctuation spectrum. The non-local component can be taken 
into account properly by considering the long range interac¬ 
tion between a pair of stripe dipoles as shown schematically 
in Figure 1. The interaction between a pair of dipoles is given 



FIG. 1. Schematic representation of the long range interaction be¬ 
tween two elementary stripe dipoles. 


by: 


SHnl 



x') 

\x — x'\ a 


(3) 


where SAi and 4.42 are the corresponding areas (see Figure 
1). If A is the modulation length of the stripe pattern, in the 
limit \x — x'\ 3> A a multipolar expansion of the interaction 


(3) leads to (see Supplemental Material): 


U n i = 7 t 7 / drx / d 2 x'VL(\x - x'\) 


e(x) ■ e(x') 


- (« + 2 ) 


x — x '|“+ 2 

e(x) ■ (x — x')e(x') ■ (x — x')' 


_ rpt I Q!+4 


\X — X‘ 


(4) 


In this expression 7 = aJP 2 and P = j J A dxx<j>(x) is the 
modulus of the dipolar moment. The unit vectors e(x) give 
the orientation of the dipoles which point along the local wave 
vector of the stripe pattern and f 2 (x) is a short range cutoff. 
Here we have neglected fluctuations in the modulation length 
and accordingly the elastic coefficient 7 is evaluated in its 
mean field value (see Supplemental Material for a discussion 
on relevant fluctuations). As we can see from the expression 
obtained, the long range repulsive interaction is responsible 
for a generalized dipolar contribution to the total energy. 

Considering again small fluctuations in the direction of the 
wave vector kg, we can write e(x) ~ ^ +Vu(x), which leads 
(considering that kg points in x direction) to: 


A Hni = ^7 [ d 2 x [ d 2 x'n(\x-x'\) 


d y u(x)dy'U(x') 


v_ _ 

_ ^/la+2 


\x — X 


- (a + 2 ) 


(y - y') 2 d y u(x)d V 'u(x') 


(5) 


Thus, the effective Hamiltonian for the displacement field 
u(x) results: 


A U=\ J ~^^{lxk 2 x + + ^nik 01 2 k y 4 )u{k)u{-k) 

(6) 

where 7 n i = 7«(2 + a) 2 C(a + 4) and the function C(a) = 
2 2 -“r(^p)/r(f) with r(x) being the Gamma function. 
From the previous considerations we are now in a position 
to analyze the stability of positional and orientational order of 
the stripe structures when the long range interactions are taken 
into account. 

Positional order: From the effective Hamiltonian ( 6 ) we 
can see that for a > 2 the 7 y k 4 term dominates over 
7 nik a ~ 2 k y 4 in the long wavelength limit, i.e. for suffi¬ 
ciently short range interaction no positional order is pos¬ 
sible. If a < 2 the term 7 n ik a ~ 2 k y 4 dominates over 
7 yk 4 , but even in this longe range interacting regime it is 
easy to check that the average square fluctuations ( u 2 ) = 

ksT J (~f x kl + 7 y k 4 + 7 nik a ~ 2 k y 4 ) 1 diverge with 
some power of the system size for any a > 0. 

Orientational order. It is well known that in systems with 
short range interactions orientational order can be weakened 
by the presence of topological defects [22, 25], The typical 
situation in two dimensional systems with continuous sym¬ 
metries is that only quasi-long range order is possible when 
interactions are of sufficiently short range [26]. Nevertheless 
it is commonly argued that even in systems with long range 
interactions, like Coulomb or dipolar interactions, shielding 
effects make the effective interactions short ranged. Here we 
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revisit this question, considering explicitly the effects of the 
range of the interactions and show that, although the shielding 
occurs, the effective interactions are still capable of stabiliz¬ 
ing a long-range-ordered nematic phase in two dimensions for 
long enough interaction range. 

At low temperatures the stripe structure can be thought of 
as composed by a mosaic of domains of average size £„ cor¬ 
responding to the correlation length of the displacement field 
u(x). The orientation of each domain is a natural order pa¬ 
rameter which can be described by a unit vector n. This vec¬ 
tor represents the mean orientation of the elementary dipoles 
inside a domain and consequently it is defined in terms of the 
unit vectors e(x) previously defined in (4) as : 


n{x) 


JAA u d2xg W 

Iaa u d2x e( x) 


(7) 


where A A u is the area of the domain, and it is over this area 
that a coarse graining process is made. Proceeding as in the 
analysis of positional order, we can separate the contribution 
to the orientational energy into two parts, a local part com¬ 
ing from interactions between nearby domains and a non-local 
one due to interactions between far apart domains. In the long 
wavelength limit, at the scale of the correlation length the 
effective interaction between nearby domains will be of the 
form: 


A n 



I 


d 2 x (X79) 2 (x) 


( 8 ) 


where 6(x) is the angle between two neighboring domains 
pointing along directions n and n'. The elastic coefficient 7 
can be estimated to be Ja 2 P 2 /( 4£“ -2 ). To continue with our 
analysis we realize that over length of order deviations of 
the local directors e{x) are small. This means after a coarse 
graining process, the interactions between far apart well po¬ 
larized domains ( of typical size £„ x ) has the same form 
of Eq. (4): 


dioni = 77 / d 2 x / d 2 X VL(\x - 0? \) 


n{x) ■ n( x') 
\x - x'\ a + 2 


— (ct + 2) 


n(x) ■ (x — x’)n{x') ■ (x — x') 


(9) 


as a consequence of the principle of superposition. Then, 
H 0 = Hoi + Honi is the complete orientational effective 
Hamiltonian. This is one of the main results of our work 
(see Supplemental Material). Note that usually the effective 
orientational energy is taken to be composed only by the lo¬ 
cal part, corresponding to smooth variations in the mean di¬ 
rections of neighboring striped domains. We will see in the 
sequel that the presence of the second (non-local) term can 
potentially change the universality class of the orientational 
order in the system. A renormalization group study of the ori¬ 
entational effective Hamiltonian H a has been done before in 
Ref. 27 for the case a = 1, which corresponds to a dipolar XY 
model. In that reference, the authors were able to renormal¬ 
ize the model and, importantly, they showed that the universal 


properties are not changed by the presence of the anisotropic 
part of the interaction. Furthermore, they showed that a whole 
family of models with isotropic long range interactions of the 
form fr | k\ a S(k)S(—k) behave in qualitatively the same way 
as the dipolar XY model as long as the range cr < 2. Once the 
mapping between these models and ours is established then 
the critical properties of the stripe forming systems are known. 
In fact, the Fourier transform of the isotropic term in Eq. (9) 
is proportional to Jr \k\ a n(k)n(—k) [ 8 ], Then, one immedi¬ 
ately see that for a > 2 the leading term in 77,„ is quadratic 
in k. In this case the low temperature physics of the system is 
that of the two dimensional short range XY model, i.e. there 
is a phase transition of the Kosterlitz-Thouless type at a crit¬ 
ical temperature Tkt- In a system with dipolar interactions 
a = 3 and then we expect it to have an isotropic-nematic 
phase transition of the KT type, as anticipated in previous 
works based on analysis of fluctuations of the local part of 
the effective Hamiltonian [23, 28]. In this case nematic order 
is quasi-long-range with algebraically decaying correlations. 
However, when a < 2 the physics changes according to the 
results of Ref. 27. Now, the non-local part in H 0 is rele¬ 
vant and rules the low temperature phase transition. In fact, 
the long range nature of the interactions in this sector are able 
to stabilize a nematic phase with truly long range order be¬ 
low a critical temperature T c . It is possible to show, in the 
framework of renormalization group equations, that the criti¬ 
cal properties of the systems for a < 2 show some peculiar 
characteristics, for example [27]: 

• in the critical region, the correlation length diverges 
exponentially at T c , from both sides, as £ 0 oc 

exp ( , = I, reminiscent of the KT transition be- 

Vvi T c- T i / 

havior. 

• For T < T c in the critical region, the average dipolar 
moment behaves as M oc S, 0 {T)~^ 2 ~ a ^ 2 , showing the 
existence of long range order when a < 2. 

• The orientational susceptibility diverges as Xo oc 
£o(T)“ in the critical region. 

This kind of behavior should be observable, e.g. in systems 
with long range Coulomb interactions for which a = 1. This 
case maps onto the dipolar XY model analyzed in 27 and 
the results may be relevant to understand the phase behav¬ 
ior of two dimensional electron systems. In the next section 
we show results from computer simulations of systems with 
a — 1 (Coulomb) and a = 3 (dipolar) which give support to 
the different scenarios in both systems as described before. 

Simulation results: We performed Langevin simulations of 
the Hamiltonian (1). The relaxational (overdamped) Langevin 
dynamics of the density <f>(x) is defined in reciprocal space by: 

t^(M) = -Aoc{k)<j){k,t) - w[0 3 ]f(M) +»?(M) (10) 

where A a (k) is the spectrum of fluctuations, i.e. the Fourier 
transform of the quadratic part of the effective Hamiltonian 
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(1), [</ 3 ]f(/c, t) stands for the Fourier transform of </f> 3 (f, t) 
and i](k, t) represents a Gaussian white noise with correla¬ 
tions < r](k,t)rj(k',t') >= (2ir) 2 2T5(k + k')5(t — t'), 
where T is the effective temperature of the heat bath. We 
worked with two forms of A a (k), the first one _4 3 (fc) = 
0,2 {k — kf) 2 — r encodes the linear dependence of the isotropic 
dipolar interaction with k, with a 2 and r constants. The sec¬ 
ond form is .4 1 (/,') = a 2 (k 2 + 2 kQ/k— 3&q) —r, corresponding 
to the Coulomb interaction proportional to 1/fc in two dimen¬ 
sional Fourier space. The parameters were chosen such as to 
have the same values of A a (k) close to the minimum at fc 0 . To 
ensure this we have set 02 = 1 for the dipolar and a 2 = 1/3 
for the Coulomb cases. In both cases we set r = 1 and ko = 1. 

For the numerical simulations we have used an implicit 
first-order scheme for the numerical integration of (10) in the 
Fourier space, a procedure that guarantees good numerical 
stability with time step dt = 0.1, as established in previous 
works [29, 30]. In the adimensional form, the periodicity of 
the stripes are set by the lattice constant dx of a 2d square grid 
with linear size L = MN, so that k = (k x /dx , k y /dx) with 
ki = 2irrii/L and dx = 7 t/M. Within this scheme, the stripe 
length span M lattice sites and the linear system size is such 
that contains N stripes. We fixed M = 11 in order to have 
smooth domain walls. 

After an estimation of the equilibration and correlation 
times from high temperature quenches, we performed slow 
cooling experiments and found that below T ~ 0.59 (T ~ 
0.45) the dipolar (Coulomb) systems find themselves in the 
low temperature phases (with orientational order) for all sys¬ 
tem sizes. Above those temperatures the configurations are in 
a state usually called liquid of stripes, where both positional 
and orientational correlation lengths are finite. So we concen¬ 
trated on equilibrium simulations for T = 0.57 (T = 0.43) for 
system sizes ranging from (12 x ll) 2 up to (66 x ll) 2 . The 
orientational order was quantified through the local director 
field v(x) = V0(x)/|V</(a;)| by measuring Q = (cos2ff(a;)) 
and its corresponding susceptibility x„, where 0(x) is the an¬ 
gle defining the local orientation of the director field. 

The previous analysis implies that, in the limit of large sys¬ 
tem sizes, the orientational susceptibility for interactions with 
a = 1 and a = 3 should be qualitatively different for T < T c . 
In the Coulomb case, the second order nature of the phase 
transition should imply that the susceptibility must be finite 
when N -X 00 . On the other hand, for dipolar interactions 
the transition should be of the KT type, implying a monotonic 
(logarithmic) increase of \o with system size, which should 
diverge in the thermodynamic limit for all T < Tkt- Re¬ 
sults for the orientational susceptibility as a function of the 
linear system size (N = L/ 11) from simulations are shown 
in Fig. 2 for the two characteristic temperatures cited above, 
corresponding to the low temperature phase of each model. 
Although computational limitations prevent us to reach very 
large system sizes, it is clarly observed that the susceptibility 
of the Coulomb system (a = 1) first grows with N but even¬ 
tually suffers a crossover and then saturates at a fixed value 


for the largest sizes. On the other hand, the susceptibility in 
the dipolar system (a = 3) shows a power law increase with 
system size, a behavior consistent with that of a KT-like crit¬ 
ical phase. Of course, we cannot conclude that \o will not 
saturate at larger TV’s, but the different trend observed in both 
systems for equivalent parameter values is a strong indication 
that the theoretical results are indeed correct. 


256 


^ 128 


64 


12 18 24 30 42 54 66 

N 

FIG. 2. Orientational susceptibility of the Coulomb (blue squares) 
and dipolar (red dots) models as function of the linear size of the 
systems in log-log scale. The full lines are power law fits with expo¬ 
nents 0.7 (blue) and 0.6 (red). 



In summary, we have shown that two-dimensional stripe 
forming systems with isotropic competing interactions can 
be classified into two universality classes: for sufficiently 
short-range interactions a Kosterlitz-Thouless transition from 
an isotropic to a quasi-long-range orientational order phase 
takes place with the well known phenomenology of defect- 
mediated phase transitions; but, for sufficiently long-range re¬ 
pulsive interactions a second order phase transition with some 
unusual characteristics drives the system from the isotropic to 
a fully long-range orientational order phase. These results im¬ 
prove considerably the understanding of the nature of phase 
transtions in stripe forming systems and may be relevant to 
a wide variety of systems, particularly the strong correlated 
regime of two dimensional “electronic liquid-crystals” phases 
and modulated phases in ultrathin magnetic films with perpen¬ 
dicular anisotropy. 
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Supplemental information 

Nature of long range order in stripe forming systems with long range repulsive interactions 

Alejandro Mendoza-Coto, Daniel A. Stariolo and Lucas Nicolao 


Here we show the fundamental steps which lead to the effective positional and orientational Hamiltonians, and discuss the 
main approximations considered. 

Our starting point is the following coarse grained Hamiltonian for a scalar field <j)(x) containing competing interactions and a 
local potential: 

H[c/)(x)\ = 1 J d 2 x (v<j)(x)^ + 1 J d 2 x J d 2 x' <j>(x)J(x- x')(j>(x') + ^ J d 2 x V(<j>(x)), (1) 

where/? = 1/ksT and V ((f)) = —|0 2 +|</> 4 . This particular form used for the local potential is not essential for the conclusions 
of our study, we chose this particular form of V{(j)) for simplicity in the numerical simulations. 

The non-local interaction considered by us corresponds to a kernel of the form J(x) = J/\x\ a , and as a consequence the 
interacting part of the original Hamiltonian can be written in Fourier space in the form: 

n t [(t){x)\ = ^ J -^^A{k)4>{k)4>{-k), (2) 

where A(k) represents the fluctuation spectrum, which is isotropic and given by A(k) = k 2 + JC(a)k a ~ 2 , where k stands for 
the modulus of the vector k. Here we have used the well known result: 

T E = C(a)k a ~ 2 (3) 

where C(a) = 2 2- °T ( 2 ) /r (|r) and .F[...] represents the Fourier transform operator (see e.g. Ref. ? where a detailed 

discussion is present, and references therein). For a < 4 the function ,4 (/,::) has a minimum at a finite value of wave vector k m ’ . 

It is known that in the class of systems described by the effective Hamiltonian (1) the competition between the short range 
attractive interaction and the long range repulsive interaction gives rise to a modulation in the order parameter in the form of 
stripes, at least at zero external field. We now proceed to analyze local fluctuations of a typical stripe pattern. 

Locally, one can assume without loss of generality that the modulation profile extends in the “x" direction. Then, the scalar 
order parameter can be written in the general form <j>(x) = X^ n >i cos{nkox), where the amplitudes </>„ and ko are parameters 
which minimze the free energy. In particular, note that ko will be in general a temperature dependent parameter proportional 
to the inverse modulation length of the stripe solution. We can now study the energy cost of perturbing such a solution by a 
fluctuation u(x) in the phase of the modulation. The modulation profile changes to: 

</>(£) = E <j> n cos(nkoX + nkou{x)). (4) 

n> 1 

In the small fluctuations limit the Fourier transform of 0 is given by: 

4>(k) = (6(k-nk 0 ) + 6(k + nk 0 j^ - E (u{k-nk 0 ) - u(k + nk 0 )J , (5) 

n> 1 n> 1 

where ko = koe x . Substituting this form of <f> back into equation (2) and neglecting zero average terms we get for the leading 
contribution: 

Hi = J E A ( nk °) 0™ + \ J (^2 u{k)u(-k), ( 6 ) 

where 

12 2 

B(k) = ^ nTl ~° (27t) 4 — nko) + A(k + nk 0 )j , 

n> 1 


( 7 ) 
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and V is the system size. In the long wavelength limit (k —> 0) the dominant contribution to the spectrum function in (7) is 
around ±rifco, where ko satisfies: 

< 8 > 

An expansion of (7) in the long wavelength limit leads to: 

B(k) = y x kl +'Yyk*. (9) 

Equations (6) to (9) imply that the local effective Hamiltonian 

Ui = H 0 i + * J ( 7 xk 2 x + 7 y ky) u(k)u(-k) ( 10 ) 

is of a more general nature than usually assumed. In particular, it is not restricted to a single mode approximation which should 
limit its validity only to a neighborhood of the critical point. This result also implies that although inclusion of higher order 
modes enhances the values of the elastic coefficients 7 X and j y , this growing stiffness will not be able to drive the system to a 
state with long range or quasi long range positional order. 

To proceed further, it is important to note that the previous expansion of B{k) in the long wave limit actually represents small 
deviations in the local wave vector of the modulation (fco + koVu(x)) from its optimum value. This means that regions far away 
from k 0 were neglected in the original spectrum A(k) in equation (2). In particular, energy contributions coming from the long 
range tail of the repulsive interaction were left out from the previous analysis. 

We go on to the calculation of the explicit interaction between far apart domains (basic cells), considering the basic cell as a 
region containing a single interface with average length equal to the modulation length A, as depicted in Figure 1. 



FIG. 1. Schematic representation of the long range interaction between two elementary stripe dipoles. 
The interaction between a pair of dipolar cells is given by: 
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In the limit of interest, \x — x'\ A, the integral can be calculated to any order by a multipolar expansion. Keeping terms up to 
second order we find: 
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where x 2 i represents x — x!, Qi = (1 /SAj) f gA d 2 x <p(x) is a generalized “charge” and Pi = (1 /SAi) f SA d 2 x x 4>(x) is a 

generalized “dipole moment”. Then, up to second order, the interaction between two basic cells is composed by a term like a 
charge-charge interaction, two terms of type dipole-charge interaction and a generalized dipole-dipole type interaction. 

The total contribution of this interaction energy in the continuum limit is: 
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(13) 


where f2(|a?2i|) is a short distance cutoff. Now, translational invariance can be invoked to conclude that ( Q(x 2 )P{xi )) = 
(Q(xi)P{x 2 )), where (• • •) means a thermal average. This implies that the “dipole-charge” interactions do not contribute to 
the free energy. In other words, charge and dipolar degrees of freedom are decoupled, and because of this the “charge-charge” 
interactions are irrelevant for the orientational order. Thus, the effective Hamiltonian with the degrees of freedom responsible 
for orientational order reduces to: 
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which is equation (4) in the letter, where e(x) = P(x)/P(x) represents the unit vector pointing in the direction of the dipolar 
moment of a given cell. The directors can be written in terms of the displacement fields u(x) as e(x) sa ^ + Vu(i) and then, 
to the leading order in u(x) Eq. (14) can be written as: 
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Up to this point, the dipolar moments in the orientational effective Hamiltonian can fluctuate both in orientation and in 
magnitude. It is worth to note that in this kind of systems a disordered stripe structure or “liquid of stripes” develops at energies 
or temperatures larger than the establishment of orientational or translational order. After phase separation has been settled 
in the form of a liquid of stripes, the stiffness of the stripes in the longitudinal direction greatly exceeds the stiffness in the 
transversal direction. Actually, this is a reflection of the existence of Goldstone soft modes in the transversal direction, evident 
in the isotropic nature of the spectrum of fluctuations at wave vector k m in the stripe liquid state, which plays a major role in 
the kind of order at low temperature in these systems. The possible order-disorder transition and phase separation in these stripe 
systems are in different energy scales. It is around the phase separation transition point that fluctuations in the width of the 
stripes (magnitude of generalized dipolar moments) are important. At lower temperatures, when the liquid stripe structure has 
already been settled, longitudinal fluctuations are very small. Then, in order to study the nature of orientational order, it seems 
justified to consider |P(x)| as being constant. In this way we arrive at the final form of the orientational effective Hamiltonian 
(eq. (5) in the letter): 
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where 7 = aJP 2 . It is instructive to analyse the non-local effective Hamiltonian in Fourier space. It is known that 
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where C(a) = 2 2 °T (Pf^) /T (|) and k = yJk 2 + k y . Then it is straightforward to conclude that T = C(a + 2 )k a , 

d 2 {k 2+a ) 


and 


T 


r a+4 


= —C(a + 4)- 


dk 2 y 


= —C{a + 4)(2 + a)k a - C(a + 4)(2 + a)afc a “ 2 A: 2 . 


(17) 


With these results the energy (15) can be written in Fourier space as: 
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Noting that C(a + 2) + C(a + 4)(2 + a) 2 = 0, the previous result simplifies to: 

A K n i = | J -^21 k a - 2 k*u(k)u(-k) (19) 

where l = C(a + 4) (2 + a) 2 a. Finally, summing the local contribution of eq. (10) with the non local one of eq. (19) we arrive 
at the effective Hamiltonian for the positional degrees of freedom (eq. (6) of the letter): 

A U = \j ~^^(^ k l + ^yK + ^ik a ~ 2 ky 4 )u(k)u(-k) (20) 

where 7 ni = 701(2 + a) 2 C(a + 4). 

This is the relevant effective Hamiltonian which determines the nature of positional order in the stripe system. The next step 
is to obtain an effective Hamiltonian for the orientational degrees of freedom. Follwing the previous method of analysis, one can 
argue that the effective orientational Hamiltonian will be composed by two contributions. The local one, that have been obtained 
in previous works, can be seen as the result of a coarse graining process over lengths of the order of the positional correlation 
length Beyond this lengthscale the director field of the stripes will show deviations which can be modelled, in the small 
deviation regime, by an effective orientational local Hamiltonian of the form: 


A -Hoi = ^ J d 2 x(V9) 2 {x) 


( 21 ) 


where 9(x) is the angle between two neighboring domains pointing along directions n and n'. Note that these stripe domains 
have linear size of order £ u , much bigger than the size of elementary dipoles considered previously for the study of positional 
order. The elastic coefficient 7 ^ can be estimated to be Ja 2 P 2 /( 4£“ -2 ). On the other hand, the interaction between domains 
separated by distances much larger than £„ can be taken into account following formally the same steps as done for the non local 
interaction between positional degrees of freedom. The effective interaction between domains will take the same form as (14): 
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where 7 = aJP 2 . Thus, A? i 0 = A Hoi + A H 0 ni is the complete orientational effective Hamiltonian. To study the stability of 
an ordered orientational phase against small deviations of the local directors n(x) a spin wave analysis can be done. In the spin 
wave limit the orientational Hamiltonian reduces to: 
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From Eq. (19) the Fourier transform is easily found to be: 

n° = \J (g, (7 oik 2 + 7 nik a - 2 k y 2 ) 9(k)9(-k) 


(24) 


with 7 ni = 7 a (2 + a) 2 C(a + 4). 



